Prom one to N Cooper pairs, step by step 
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We extend the one-pair Cooper configuration towards Bardeen-Cooper-Schrieffer (BCS) model of 
superconductivity by adding one-by-one electron pairs to an energy layer where a small attraction 
acts. To do it, we solve Richardson's equations analytically in the dilute limit of pairs on the one- 
Cooper pair scale. We find, through only keeping the first order term in this expansion, that the N 
correlated pair energy reads as the energy of N isolated pairs within a N(N — 1) correction induced 
by the Pauli exclusion principle which tends to decrease the average pair binding energy when the 
pair number increases. Quite remarkably, extension of this first-order result to the dense regime 
gives the BCS condensation energy exactly. This leads us to suggest a different understanding of 
the BCS condensation energy with a pair number equal to the number of pairs feeling the potential 
and an average pair binding energy reduced by Pauli blocking to half the single Cooper pair energy 
- instead of the more standard but far larger superconducting gap. 
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I. INTRODUCTION 

The possibility for an electron pair with up and down 
spins to form a bound state under a very small attractive 
potential has been pointed out more than 50 years ago, 
first by Frohlichi and then revealed by Cooper— who con- 
sidered adding a pair of electrons with opposite spins and 
zero total momentum, to a full Fermi sea in an energy 
layer where a weak attraction acts (see Fig. 1(a)). This 
attraction leads to the bound state for the two added 
electrons known as a " Cooper pair" . In a milestone pa- 
per of the following year, Bardeen, Cooper and Schri- 
effer (BCS)^ have succeeded to extend this idea to N 
pairs. This was a crucial step toward the understanding 
of superconductivity. The configuration they considered 
is shown in Fig. 1(b): the attractive potential now acts 
in a layer extending on both sides of thenoninteracting 
electrons Fermi level and pairs are formed from electrons 
located in this layer. This BCS configuration can be for- 
mally achieved by extending Cooper's model, through 
adding more and more pairs to the "frozen" Fermi sea 
until this layer becomes half-filled. The major difficulty 
to go from 1 to N pairs, step by step, comes from the 
proper handling of the Pauli exclusion principle between 
a fixed number of indistinguishable electrons. A smart 
way to overcome this difficulty is to turn from the canon- 
ical ensemble in which N is fixed, to the grand canonical 
ensemble in which the mean value of N is adjusted in 
the end to fit the particle number. These two procedures 
are known to exactly give the same results in the large 
TV limit - which is the physically relevant limit for the 
problem at hand. For a discussion of the applicability 
of the canonical ensemble to superconductivity, see, e.g., 
Refe4£. 

A few years later, R.W. Richardson^"— has proposed a 
way to stay in the canonical ensemble^. Through a el- 
egant procedure, he succeeded to analytically transform 




Figure 1: Fig. 1(a). Frozen Fermi sea \Fo) (grey area) in 
the one-Cooper pair problem and also in the problem we here 
address, the number of pairs added to \Fo) increasing, one 
by one. A small attractive potential extends over an energy 
layer SI above the Fermi sea \Fo). The two dots schematically 
show the two electrons with opposite spins and momenta, of 
the one-Cooper pair problem. 



the Schrddinger equation for N pairs into a coupled set 
of N algebraic equations. These equations, which are 
commonly called Richardson's equations (or more prop- 
erly Richardson-Gaudin's 10 equations), have been solved 
analytically in a few cases only: In Refj£, R. W. Richard- 
son solved them for N pairs in one and two multiply- 
degenerate levels. In another work4, he used this pro- 
cedure in the limit of an infinite number of pairs within 
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Figure 2: Fig. 1(b). Normal Fermi sea \F) (grey area) in 
the usual BCS model. The potential extends symmetrically 
(+u>c, -w c ) on both sides of this Fermi sea, so that its total 
extension is 2uj c = £1. The Fermi sea \F) of Fig. 1(b) contains 
TV = pouic — pofl/2 more electrons than the frozen Fermi sea 
\Fo) which is here shown as a dashed circle, the density of 
states po being taken as constant over fl. 



the half-filled configuration and succeeded to rederive the 
BCS grand canonical result at the leading order in the in- 
verse particle number. Nowadays, Richardson's approach 
is widely usedli for the analysis of the superconduct- 
ing properties of small granules with nanometer sizes. 
Such granules accommodate countable number of corre- 
lated electron pairs for which Richardson's equations can 
be solved numerically. Conceptual difficulties in solv- 
ing Richardson's equations can be attributed to the fact 
that this procedure essentially corresponds to reduce the 
quantum many-body problem to a classical one for "par- 
ticles" interacting via Coulomb-like forces. It however is 
well-known that the analytical treatment of a many-body 
problem, even with classical particles, stays difficult. In 
addition, a severe singularity problem arises when solv- 
ing these equations numerically^, due to possible cance- 
lations in the denominators. This is why the resolution 
of Richardson's equations in the general case, stayed an 
open problem up to now. 

The purpose of the present work is to extend the orig- 
inal Cooper's work2, from 1 to AT pairs within the canon- 
ical ensemble by solving Richardson's equations analyt- 
ically. We here calculate the ground state energy of N 
pairs when /V increases step by step starting from N = 1. 
We have recently derived the energy of just two pairs by 
various procedures^ to point out the trend induced by 
Pauli blocking on Cooper pairs in the simplest configu- 



rations. We have also given a sketch of the derivation 
for general TV in RefJ^ where the analytical expression 
of the energy of N Cooper pairs was for the first time 
reported. Although, at the present time, it seems diffi- 
cult to perform such a pair increase experimentally, this 
study can at least be seen as a gedanken experiment to 
reveal a possible connection between two famous prob- 
lems, namely, Cooper's model and the BCS model for 
superconductivity, in order to more deeply understand 
the role of the Pauli exclusion principle in Cooper-paired 
states. 

From a careful analysis of the structure of Richard- 
son's equations for a number of pairs going from 2 to 7, 
we have succeeded to see how these equations can be ma- 
nipulated to get the A^-pair energy En analytically for 
N = 2n and 2n + 1, the structure of the solution depend- 
ing on the parity of N. The analytical solution we pro- 
pose relies on an expansion in a dimensionless parameter 
proportional to the inverse of the sample volume, through 
a procedure similar to the one we have recently used to 
get the energy of two pairsi£. This parameter actually 
is the inverse of the number of pairs N c from which sin- 
gle Cooper pairs would start to overlap. This limitates 
the validity of the expansion we perform to the dilute 
regime, i.e., A" small compared to N c since we moreover 
only keep the lowest order term in this expansion. This 
restriction a priori prevents the extension of our result to 
the dense BCS regime in which half the potential layer is 
filled by electrons. It however lets N to be arbitrary large 
in the thermodynamical limit because only N/N c mat- 
ters. In spite of this limitation, it is quite remarkable to 
note that the analytical result we obtain, extended out- 
side its mathematical validity range, produces the BCS 
condensation energy exactly, i.e., with all its numerical 
prefactors, even when the potential does not extend sym- 
metrically on both sides of the Fermi level. This leads us 
to think that our solution, although derived in the dilute 
limit, can be considered as an important step towards 
mastering Richardson's approach in the general context 
of superconductivity. 

The expression we find for the ground state energy of 
N pairs casts as 

E N = NE 1 + N{ - N 2 l) E lnt (1) 

E\ is the energy of one isolated pair, as originally found 
by Cooper, while Ei nt can be seen as the interaction en- 
ergy between two pairs. We wish to stress that this in- 
teraction energy can only come from the Pauli exclusion 
principle between the electrons from which these pairs are 
constructed, because of the very peculiar form for the in- 
teraction potential taken in the BCS theory of supercon- 
ductivity which prevents any direct interaction between 
pairs. Such a N(N—l) term may appear as natural in the 
dilute limit. However, no other terms in N(N — 1)(N — 2) 
and higher seem to exist in the dense BCS regime, as 
the agreement of our dilute limit result with the dense 
BCS one tends to indicate. This agreement, highly un- 
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expected at first, brings some new light to the concept of 
" Cooper pair" . In particular, it leads us to reconsider the 
standard understanding of the BCS result for the ground 
state energy of the superconducting state in terms of the 
excitation gap. 

(i) The BCS result for the condensation energy (see, 
e.g., Refjiir— ) is usually understood through the con- 
cept of "virtual pairs" , as introduced by Schrieffer. These 
pairs are made of electrons with opposite spins and zero 
total momentum, excited above the Fermi level of non- 
interacting "normal" electrons. From the energy scale 
A over which the pair density changes in the vicinity of 
this Fermi level, one deduces a "virtual pair" number of 
the order of the number of pairs within a thin layer of 
width A. Such a number is far smaller than the total 
number N of electrons feeling the attractive potential 
and used in the BCS wave function ansatz. In order to 
recover the BCS expression of the ground state energy 
with such a small "virtual pair" number, it is necessary 
to give to each of these pairs an energy of the order of 
A. It however is clear that this "virtual pair" picture 
is highly qualitative: these pairs are not introduced ab 
initio in the BCS theory and their number is not rigor- 
ously controlled, so that the binding energy of a virtual 
pair is not a precisely defined quantity. These virtual 
pairs essentially are a qualitative concept which provides 
a physical picture for the BCS ground state energy- ^ 14 ' 15 . 
It actually is possible to introduce these excited pairs in 
a well defined manner, as done by Leggett 1 ^ using de- 
struction and creation operators acting not on the frozen 
Fermi sea, but on the normal electron Fermi sea. How- 
ever, this procedure does not yet generate a well defined 
pair binding energy. 

(ii) Our alternative understanding of the BCS result for 
the superconducting ground state energy is based on the 
fact that correlations exist between all electrons which 
feel the attraction, so that all these electrons must form 
pairs - not just a fraction as with the virtual pair un- 
derstanding described above. These pairs actually corre- 
spond to what Schrieffer— calls " superfluid pairs" . The 
pair number N being then much larger than within the 
usual BCS understanding, the average binding energy of 
these pairs ends by being much smaller than the energy 
gap A. We here show that it only is a fraction of the 
energy of one isolated Cooper pair. More precisely, this 
average binding energy decreases linearly when the num- 
ber of pairs increases. This decrease entirely comes from 
the Pauli exclusion principle that restricts the number of 
available excited states when constructing paired states: 
with less and less empty states enjoying the attractive 
potential, the correlated pairs gain a smaller and smaller 
average binding energy, as if each of the empty states 
were bringing a constant elementary contribution to this 
binding energy. Our present work actually shows that, 
in spite of the widely spread idea that the isolated pair 
picture has little meaning in the dense regime 1 ^, it is pos- 
sible to recover the BCS ground state energy in an ex- 
act way, by simply using the concept of isolated Cooper 



pair but by properly including Pauli blocking between 
the electrons of these pairs. 

We also would like to mention that the result for the 
ground state energy we found, is in agreement with the 
solution of Richardson's equations in the large-sample 
limit obtained earlier—, except that we here consider 
arbitrary fillings of the potential layer, in contrast to 
Refii which focuses on the half-filling configuration only. 
Changing the layer filling enables us to explicitly estab- 
lish the link between the Cooper model and the BCS 
picture, as they are presented in textbooks. Moreover, 
our procedure allows us to reveal an important under- 
lying physics which is somehow masked when consider- 
ing only one particular configuration, namely, half filling 
this configuration being, actually, quite specific. In ad- 
dition, the analytical approach of Refi 7 -, as well as more 
recent studie d 18 ' 19 based on the same approach, assumes 
that Richardson's energy-like quantities in the ground 
state are organized into an arc in a complex plane. This 
very strong assumption is based on some numerics on 
Richardson's equations 7 . However, in our opinion, this 
assumption is not fully controllable, so that it is not 
guaranteed that it does not break in certain region of 
the phase diagram (for instance, in the delicate BEC- 
BCS crossover region). Therefore, alternative methods 
for solving Richardson equations, as the one we here pro- 
posed, are highly desirable. 

The present paper is organized as follows. 

In Section II, we settle the problem. We outline 
Richardson's procedure to transform the Schrodinger 
equation for 2 Cooper pairs into two nonlinear equations 
and we report how this procedure extends to N pairs. We 
then explain how an analytical solution of these Richard- 
son's equations can be obtained, thanks to the identifi- 
cation of a small dimensionless parameter. In practice, 
this is done through the expansion of sums similar to the 
one which appears in the one-Cooper pair problem. 

In Section III, we use this expansion to explicitly solve 
Richardson's equations analytically, step by step, from 2 
to 7 pairs. We determine the N energy-like quantities 
which are solution of these equations, at the lowest order 
in the dimensionless parameter we have previously iden- 
tified. A careful analysis of these solutions allows us to 
construct a general form for the energy of an arbitrary 
number of pairs. The structure of the equations for In 
and In + 1 being different, the solution for 6 and 7 pairs 
corresponds to n = 3 which is somewhat generic. In con- 
trast, n = 2 for 4 and 5 pairs is not because n and n! 
are then equal. This explains why we have had to go up 
to (6, 7) to fully control how the general solution devel- 
ops. In this section, we also discuss the validity of the 
procedure we propose to solve these equations. 

In Section IV, we reconsider the condensation energy 
obtained by Bardeen, Cooper, Schrieffer and its usual 
understanding in terms of " Cooper pairs" , in the light 
of the result we have obtained through the Richardson's 
canonical ensemble procedure. 

We conclude in Section V. 
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II. EXACT FORM OF THE iV-PAIR ENERGY 
A. The problem 

We consider a "frozen Fermi sea" |F ) made of Np free 
electron pairs with up and down spins, their energy ex- 
tending from zero to ef q = /2m, and we start adding 
zo-momentum electron pairs to this Fermi sea in order 
to extend Cooper's model towards the traditional BCS 
configuration. 

Electrons outside |F ) feel an attractive potential V 
that we take as being the standard potential used in the 
BCS theory of superconductivity, without questioning its 
oversimplifications since we want to position the present 
work within the very first version of BCS theory. This 
potential reads as 



V = ^ ^k'k4't a -k'+«-Ha k t 

k,k' 

= E^'^k (2) 



k,k' 



T4;'k is further on approximated by a contact separable 
potential 



(3) 



with V being a positive constant while »k = 1 inside 
the potential layer, i.e., for sf < £k < £_f + ^ an d 
otherwise. Thus, this potential is only felt by electrons 
with opposite spins and opposite momenta, outside the 
Fermi sea |F ), up to a cut off that we have chosen to 
call Sf The system Hamiltonian then reads as H = 
H o + V where 



J2 ek ( a i 



kt°kt + a ki a H 



(4) 



We add N = 1, 2, 3... pairs of electrons with up and 
down spins and zero total momentum to this Fermi sea 
\F ) and we look for the exact ground state energy of 
these N pairs resulting from both the potential V and 
the Pauli exclusion principle. 



B. One Cooper pair 

The N = 1 case is the standard problem studied by 
Cooper. By writing the wave function as 



l^) =J2^G(k)at t al u \F Q ) 



(5) 



where \Fq) is the previously defined "frozen Fermi 

sea", we find that fulfills the Schrodinger equation 

(H — Ei) \tpi) — 0, provided that the G(k)'s are such 
that 



G(k) (2e k -E^-VY, ^pG(p) = 



(6) 



since u>k = w^.- If we now extract G(k) from the above 
equation, multiply it by Wk and sum over k, we end with 
the textbook equation for the one-Cooper pair energy, 
namely 



2e k - E 1 



(7) 



To calculate this one-pair energy explicitly, we turn to 
the continuous limit. If the density of states p(e) does not 
change too much over the potential extension ef < £k < 
e F „ + SI, we get 



PoV 
2 

(v/2)ln 



2e, 



2de 
2e-E 

, + 2n - e 



2ep — E 



(8) 



where po is the average density of states in the potential 
layer, while v = pqV is the small dimensionless param- 
eter which characterizes the weak potential felt by elec- 
trons above £F -When inserted into Eq. (7), this gives 
the binding energy of one "Cooper pair" as 



Ei — 2ef q — e c 
e c = 2fly c 



(9) 



where y c is an extremely small positive constant related 
to v through 



= ,-2/ 



Vc = e 



-2/v 



(10) 



in the weak coupling limit, i.e., v <^ 1. This shows that 
the binding energy of one Cooper pair depends in a singu- 
lar nonpertubative way on the small attractive potential 
V. This also shows, more importantly for the problem 
we here address, that this binding energy depends lin- 
early on the potential extension 0, i.e., on the number of 
empty electron states feeling the attractive potential. To 
evidence this understanding, we can rewrite e c as 



(pofi) 



(/soft) 



(11) 



where Nq — pqQ is the number of empty pair states for 
electrons with up or down spin in the energy layer Q 
while ey given by 



jle-2/poV 



(12) 



can be understood as the contribution of each of these 
empty pair states to the pair binding energy. This remark 
is going to be quite important for the understanding of 
what really happens when the number of added pairs in- 
creases to possibly fill half the potential layer - as in the 
usual BCS configuration - or even the whole layer, which 
is an interesting mathematical limit. Actually, a binding 
energy proportional to the number of empty states Nq, 
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as in Eq. (11), already makes us guessing that the av- 
erage binding energy in the iV-pair configuration should 
decrease linearly when the number of pairs increases, due 
to the decrease of empty states available in the linear 
combination responsible for the single-pair binding en- 
ergy. With this in mind, we now study 2 pairs, then 3 
pairs and so on, in order to understand how the energy 
changes with N. 



C. Two Cooper pairs and the Richardson's 
procedure 

Let us first rederive Richardson's equations for the 
reader not fully aware of them since the present work 
heavily relies on these equations. The derivation we here 
propose is somewhat different from the original one. 

The very peculiar form of the potential V given in Eq. 
(2) leads us to look for the ground state of two pairs as 

|V>2> = J2^^G(^iM)PiA 2 \F ) (13) 

with G(ki,k2) = G(k2,ki), as always possible to en- 
force due to bosonic symmetry in fermion pair exchange, 

namely, < = 4 A ■ 

We then note that for two free pairs, we have 



This makes G(pi,p2) undefined for pi = p2, which is 
unimportant since the state J3^ \F ) reduces to zero 
for ki = k2 anyway. In contrast, G(ki,k2) for ki ^ k2 
must fulfill the integral equation i ;l (k 1 ,k 2 ) = 0. The 
really unpleasant part of this integral equation is the 
coupling between (ki, k 2 ) through (2eki + 2£k 2 — E 2 ). 
To decouple k x from k 2 , Richardson suggested splits E 2 
as Ri + R 2 . In view of the structure of the equation 
F(ki,k2) = 0, we are led to find a symmetrical func- 
tion of (ki, k 2 ) which, multiplied by (26^ + 2ek 2 ~ E 2 ), 
splits into a ki function plus a k2 function. To do it, we 
first note that this unpleasant factor also appears in the 
identity 



2£ki + 2£k 2 



2 £kl -i?i 2e k2 -^2 (2e kl - i?i) (2e k2 - R 2 ) 



(19) 



and a similar identity with (R 1 , R 2 ) exchanged. Since 
we look for G(ki,k2) as a symmetrical function of (ki, 
k 2 ), we are naturally led to consider 



G(ki,k 2 ) = 



1 



(2eki — Ri) (2ek 2 — R2) 



(2e kl - R2) (2e k2 - Ri) 



(20) 



HnfllA. = ^ + 2£k 2 ) PlPl \F Q ) (14) 



t at 



when multiplied by (2ek x + 2ek 
resulting quantity splits 

plus a similar term 



E 2 ), 
into 



V/^A \Fo) 



Indeed 
the 

(2e kl -R 1 )~ 1 + {2e kl - Ri 

with ki replaced by k 2 . We then find that this G(k 1; k 2 ) 
function fulfills Eq. (16), provided that 



- V(l <5 kl k 2 ) (wkJI E^p l F °> ( 15 ) w ^ ( 2£k \ Rl + 2£k ^ ) = Vw ^ E ^G(k,kx) 

P 1 k#ki 



the prefactor (1 — <5kik 2 ) m the RHS of the above equation 
being necessary to enforce the cancellation of this RHS 
for ki = k2, the state /3 k 2 \Fq) then reducing to zero. As 
G(ki,k2) = G(k2,ki), this leads to 

= (h e 2 ) iv> 2 > = E F ^iM)vlA. i^o> (i6) 

where we have set 

F(ki , k 2 ) = w kl w^ 2 [(2e kl + 2e k2 - E 2 ) G(ki , k 2 ) 

- ^ E wkG ( k ' ki ) - v E wkG ( k ' k2 )i 

k^ki k#k 2 

By using the scalar product of two-free pair states 

(Fo|/3 P2 /3 Pl /3 ki /3 k2 |^o) = (1--WJ 

(£ Pl ki<W 2 + 5 Pl k 2 ^ P2 kJ 

we find that \tp 2 ) is solution of the Schrodinger equation 
provided that 



(21) 

for arbitrary ki . The next step is to get rid of k x . For 
that, we first insert Eq. (20) for G into the above equa- 
tion. This leads to 



2e kl - Ri 

U>k! 



^E 



2e k - Ri 



+ 



2e kl - R2 



2ek — Ri 



2w kl 



(2e kl - Ri) (2e kl - R 2 ) 



(22) 



where the sums are now taken over all k. By noting that 
the RHS of the above equation also reads 



1 



2wk! 

Ri — R 2 \2ski — ^1 2e 



1 



?ki 



Ri 



for i?i ^= R 2 , as always possible to enforce since, up to 
now, the only requirement is to have (i?i, R 2 ) such that 
E 2 = Ri + R 2 , we then find that Eq. (22) leads to 



(l-<W 2 )F( Pl ,p 2 ) =0 



(18) 



f(Ri,R 2 



| f(R 2 ,Ri) =Q 



2£ki — ^1 2eki — R 



(23) 
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where f(Ri, R2) is given by 
f{R l ,R 2 ) = l-VY J 



2V 



2ek — i?i i?i — i?2 



(24) 



Since R\ ^ R2, the validity of Eq. (23) for arbitrary ki 
imposes f(Ri,R 2 ) = = f(R 2 ,Ri). So that we end with 



p 



+ 



2F 



2e p - i?i i?i - R 2 



w r . 



+ 



2V 



2e p — i?2 -R2 — -R 



(25) 



These two equations are known as the Richardson's 
equations for two pairs. They transform the resolution of 
the Schrodinger equation for two pairs into the resolution 
of two coupled algebraic equations for (Ri, R 2 ), a far 
simpler problem, dcfinitly. We wish to stress that the 
existence of an exact solution to the 2-pair Schrodinger 
equation is deeply linked to the very peculiar form of the 
reduced BCS interaction potential, the zero-momentum 
pair staying unbroken by this interaction. 



D. Richardson's equations for N Cooper pairs 

We now consider the Fermi sea \Fq) plus three electron 
pairs with up and down spins. A similar procedure leads 
us to split the 3-pair energy as E$ = Ri + R 2 + R3 , with 
(i?i, i?2, R3) fulfilling the three equations 



1 = V J2 
p 

p 



2V 



2V 



2e p - 


Ri 


+ 


Ri — R2 


+ 


R\ — i?3 


w p 




+ 


2V 


+ 


2V 


2e p - 


R2 


R 2 - Ri 


R2 - Rs 






+ 


2V 


+ 


2V 



2Sp — -R3 i?3 — i?i i?3 — i?2 



(26) 



And so on for a higher number of pairs. 



E. Expansion of the relevant sum 

The unpleasant parts of these equations are the sums 
over p. It is in fact possible to tackle this set of nonlinear 
equations analytically by expanding these sums close to 
the one-Cooper pair energy E\ defined in Eq. (7), since 
for R = E\, the sum reduces to 1/V. A convenient way 
to do it, is to note that 



^E 
p 

^E 



2sp - R 



R-E 1 



2e p -E 1 (2e p - E{) (2e p - R) 

OO 



(27) 



The prefactor a m in this expansion reads in terms of the 
one-pair energy E\ as 



^ (2e p -E 1 ) m+1 2m <* 



(28) 



where e c is the binding energy of one Cooper pair defined 
in Eq. (11), while I m is found to read in terms of y c 
defined in Eq. (10) as 



1 



IJc 



1 — e -2 " 1 /'" 



(29) 



This makes it very close but not exactly equal to 1. Dif- 
ference with 1 is going to be crucial in the final expression 
of the Af-pair energy. In view of Eq. (28), we are led to 
rescale the energy difference in Eq. (27) as 



R — Ei— ze c 



(30) 



This allows us to ultimately write the sums appearing in 
Richardson's equations in a dimensionless form as 



00 T 

^ 2e p - R ^ rn 



(31) 



m— 1 



this expansion being a priori valid for \z\ small compared 
to 1. 



III. ANALYTICAL SOLUTION OF 
RICHARDSON'S EQUATIONS 

We now turn to the resolution of Richardson's equa- 
tions starting from two pairs, as given in Eq. (25). 



A. Two Cooper pairs 

By using the expansion given in Eq. (31), we can trans- 
form the two Richardson's equations for 2 pairs as 







y ^ zr + _]_ 

m— 1 



TO 



7 



m — 1 



Z2 — Zl 



(32) 
(33) 



where the z^s are related to the R^s through Eq. (30), 
namely, Ri — E\ = zit c . The dimensionless parameter 7 
in these equations appears to be 



7 : 



Po£ c 
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(34) 



m— 1 



where N c = poe c is the number of states within a one- 
Cooper pair binding energy layer at the frozen Fermi sea 
level. Although small compared to the number of pairs 
Nq in the energy layer where the potential acts, N c is 
definitely very large compared to 1. Indeed, in a sample 
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volume L D . the density of states po at the Fermi level 
with momentum kp is of the order of L D mkp~ 2 , where 
D = 3 or 2 is the space dimension. This gives N c of the 
order of L D /l 2 lp~ 2 where lp — kp 1 is the distance be- 
tween Fermi sea electrons and l c , defined as e c ~ l/m c l 2 , 
scales as the distance between the two electrons forming 
one Cooper pair. This makes N c infinitely large com- 
pared to 1 in the infinite-size sample limit. Consequently 
7, defined in Eq. (34), is a small parameter for macro- 
scopic samples. 

The existence of such a small dimensionless parameter 
is the key for solving Richardson's equations analytically. 
Indeed, in order for the set of Eqs. (32,33) to have a 
solution with Zi small, as necessary for the expansion 
(31) to be valid, we must have Z\ — Z2 of the order of y^y. 
By replacing one of the two equations (32,33) by their 
sum, namely 



' « m 



(35) 



m— 1 



we find from the m — 1 term of this equation that z\ + z 2 
must be equal to zero at the y/j order. This imposes to 
look for z\ and Z2 as z\ ~ —z-i — ai^/7- When replaced 
into Eq. (32) in which we only keep the m = 1 term, we 
find ~ hai + l/2a\\ so that of ~ — 1/2 Ji. This readily 
shows that a\ is imaginary at the y/j order. 

To get a nonzero value for the sum z\ + z 2 , as necessary 
to find the energy change induced by interaction between 
Cooper pairs, i.e., to have E2 not exactly equal to 2E\ 1 
we must keep two terms in Eq. (35). This gives 



0~/ 1 (z 1 +Z 2 ) + y {zl+zl) 



so that {z\ + Z2) to lowest order in 7 reads as 



zi + 



I 2 , 2\ 1 ^ 

Z2 "-2>T (zi + Z2) " 2 7 7F 



(36) 



(37) 



If we now note that (e c j) is just (4/po), as seen from Eq. 
(34), the above equation gives the energy of two Cooper 
pairs E 2 = Ri + R2 as 



Ik 

Polf 



E1 + -4 



(38) 



Since I2/I 2 — 1 + 2y c , while 2y c po = e c /Nn, this energy 
appears in terms of physically relevant quantities as 



E 9 











PoJ 


\2e Fo - 




Pa J 



1 

2(N n - l)e v 



(39) 



As previously introduced, Nq = pofl is the number of 
pair states in the potential extension ft above the Fermi 
sea \Fo), while ey is the elementary contribution to the 
pair binding energy coming from each of the empty pair 



states feeling the potential, their number being reduced 
to (Nq — 1) in the 2-pair configuration, due to the added 
pair . 

The above equations show that the energy change from 
the energy 2E\ of two isolated Cooper pairs has a naive 
term 2/po which simply comes from the increase of the 
Fermi sea level when a second pair is added. It also has 
a more subtle term 2e c /iVo = 4ey which can be seen as 
twice the shrinkage of the one-Cooper pair binding en- 
ergy coming from the Pauli exclusion principle between 
electrons with up spin and between electrons with down 
spin, available to form Cooper pairs within the potential 
extension fi. Such a "moth-eaten effect" which here ex- 
ists for Cooper pairs, has already been found in the case 
of composite bosons like Wannier or Frenkel excitons 2 ^,. 

Let us now see how this energy change generalizes when 
the number of added pairs increases from N = 2. Since 
the analytical resolution of Richardson's equations has 
raised problems for decades, we are going to consider N = 
3,4,... pairs successively in order to, in a precise way, 
show how the analytical solution for general N develops. 
This will in particular, allow us to see that, although the 
structure of the Ri solutions is different for odd and even 
TV, the sum of R^s casts exactly in the same way. 



B. Three Cooper pairs 

We now consider the three Richardson's equations 
given in Eq. (26). We first replace one of them by their 
sum which reads as Eq. (35) with three terms instead 
of two. By keeping the m = 1 term only in this sum, 
we find zi + z 2 + Z3 ~ 0. Since, in the one hand the 
Zi — Zj differences have to be in y/j for the Eqs. (26) to 
have a solution with small, while in the other hand, 
(zi,Z2,Z3) must have some internal symmetries, we are 
led to look for a solution which reads at the lowest order 
in 7 as 



z\ = -z 2 ~ 01^7, z 3 ~ 



(40) 



When inserted into the first and third equations (26) in 
which the first term only in each sum is kept, we find that 
these two equations reduce to the same equation which 
reads as 



~ hat 



1 1 

a\ 2a\ 



(41) 



This readily gives a\ = — 3/2Ii which again makes ai 
imaginary at the lowest order in 7. If we now keep two 
terms in the sum of Eqs. (26), namely 



~ h{Zl +Z 2 + Z S ) + y { Z \ + z\ + zfj 

we find the sum of z^s at the lowest order in 7 as 

3 h 
2 1 1 2 



(42) 



(43) 
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When used into the energy of three pairs £"3 = R\ + i? 2 + 
R3, we find, since e c j = 4/ 'p , this energy as 

S3 = 3 (^ 1 + ^l) 

When compared to Eq. (38), this shows that the energy 
change per pair is just twice as large for three pairs than 
for two pairs. Let us now see how this nicely simple result 
extends when the number of pairs keeps increasing. 

C. Four and five Cooper pairs 

(i) In the case of four pairs, we have four equations 
similar to Eq. (25). We replace one of them by their 
sum. Its lowest order term gives z\ + Z2 + Z3 + Z4 — 0. 
Since the Z, — Zj differences must be of the order of y/j for 
the Zi's to be small, the expected symmetrical solution 
must read as 

zi = -z 4 ~ ai-y/7, z 2 = -z 3 ~ a 2 ^7 (45) 

When inserted into the three other equations, we find 
that two of them are identical; so that we are left with 
two equations to get ai and a 2 , namely 

1 1 1 

~ hai + + ■ + — 

ai — C12 fli + a 2 2ai 

This imposes the prefactor of 1/ai to be zero, with a sim- 
ilar result for l/a 2 . If we now add these two prefactors, 
we get 

~ h (a\ + a\) + 2 + 1 (47) 

This gives a\ + a\ — — 3/Ji. When inserted into the sum 
of Richardson's equations, with now two terms kept in 
the sum, we find z\ + z 2 + 23 + 24 ~ 37/2/^1- Since 
e c 7 = 4/po, the energy of four pairs ultimatly appears as 

£4 = 4 (e 1 + -^) (48) 

(ii) For five pairs, the sum of the five Richardson's 
equations taken at lowest order in 7, gives z\ + z 2 + 23 + 
z 4 + z 5 ~ 0. Since the differences Zi — Zj must be in 
y^y for the Zj's to be small, we expect the symmetrical 
solution to read as 

zi = -z 5 ~ ai^/7, z 2 = -Z4, ~ a 2 ^7, z 3 ~ (49) 

The remaining four equations are then identical two by 
two, so that we are again left with two equations only for 
a\ and a 2 . They read 

1111 

~ JiOi + + — + ; + — 

ai — a 2 ai ai + a 2 la\ 

( T 2 2a i 3\ 1 
= ha\ + -^-^ + - — 50 

\ a\ - a,2 2/ ai 



and similar one for a 2 . By adding the 1/ai and l/a 2 
prefactors in the above equation, we get 

~ h {a\ + 4) + 2 + 3 (51) 

instead of Eq. (47). So that we now have a\+a% = — 5/Ji. 
When inserted into the sum of Richardson's equations 
in which we now keep two terms in the sum, we find 
z\ + z 2 + z 3 + Z4 + z 5 ~ 57/2/^1! so that the energy of 
five pairs appears as 

E5 = 5 ( El + 7o^> (52) 

Although the numerical prefactors in the z$ sums do 
not seem to have any physically relevant meaning, the 
results given in Eqs. (38), (44), (48), and (52) for the 
energy of N pairs clearly indicate that the energy change 
per pair depends linearly on the number of pairs as (N — 
1), in spite of the fact that the precise structure of the 
solution of Richardson's equation is definitly different for 
N = 2n and N = 2n + 1. In order to better grasp how 
this nice analytical result shows up, let us briefly consider 
one more set of TV's, namely N = 6 and 7, since these 
correspond to n = 3 which usually is generic while n = 2 
in the case of N = (4, 5) is usually not, 2 and 2! being 
equal. 

D. Six and seven Cooper pairs 

(i) For six Cooper pairs, the symmetrical solution is 

z\ = -z 6 ~ 01^/7, z 2 = -z 5 ~ a 2v /7, z 3 = -z 4 ^ 03^7 

(53) 

Among the remaining five Richardson's equations, we 
only get three different ones for a\, a 2 , a 3 . They read 

~ ha x + + 

ai — a 2 a\ — a 3 

1 1 1 



ai + a 3 ai+ a 2 2a\ 




and similarly for a 2 and a 3 . By adding the 1/ai prefac- 
tors, we end with 

~ h {a\ + a\ + o§) + 2 + 2 + 2 + 3^ (55) 

This leads to z\ + ...+ z 6 ~ (15/2)7/2/^1! so that since 
e c 7 = 4/po, we get the expected result, namely 

£7 6 ~6(£7i + -§ > ) (56) 

(ii) For seven Cooper pairs, the symmetrical solution 

is 

zi = -z 7 ~ ai-y/7, z 2 = -z 6 ~ a 2v /7, z 3 = -z 5 ~ a 3% /7 
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z 4 ~ 



(57) 



When inserted into the five remaining Richardson's equa- 
tions, we only get three different ones. They read as 



n r 1 111 1 

~ 7iffli + + + — + ; + 



1 

a\ — a 2 a\ — C13 ' a\ ' a\ + C13 a± + a 2 2a\ 

The sum a\ + a\ + a\ ends by reading as Eq. (55) except 
for the last 1/2 factor which is now replaced by 3/2. This 
leads to z\ + ...+ z 7 ~ (21/2) 7/2/^1 • So that we find the 
result we expect, namely 

^ 7 ( Si+ ^!) (59) 



E. 2n and 2n + 1 Cooper pairs 

The above analysis of Richardson's equations for N = 
1,2,3, 7 has made us knowledgeable enough to now 
face the general N solution through a similar procedure. 

(i) When the number of pairs is even, TV = 2n, we 
expect the symmetrical solution for the Zi's at the lowest 
order in 7 to be 

z\ = -z 2n — aiyfy, z 2 = -z 2n -i — a 2 ^/j, —, 

z n = -z n +i ^ a n ^/l (60) 

By inserting these Zj's into Richardson's equations, we 
get n different equations for 01, a n . They read as 



~ hai 



1 



1 



ai — a 2 ai — a 3 



+ ... 



1 



1 



^ 77- (61) 
ai + o 3 a\ + a 2 ia\ 

When multiplied by a\ and added to similar ones for a 2 , 
a n , these equations lead to 



~ h (af + ... + a 2 n ) + 2Cl + - 



(62) 



where C\ — n(n — l)/2 is the number of ways to choose 
a couple (ai — aj) among n coefficients aj. Note that 
the above equation agrees with Eq. (47) for n — 2 and 
with Eq. (55) for n = 3. When inserted into the sum of 
Richardson's equations with now two terms kept in the 
sum, namely 



2n T 2n 



(63) 



1=1 



i=l 



we find from Eqs. (60, 62) that, at lowest order in 7, the 
sum of zi is given by 

Y,^- \Mn - l)/2 + n/2] lh/ll ^ n(n - l/2) 7 -| 
i=i 1 

(64) 



So that the sum of energy-like quantities introduced by 
Richardson appears to read as 



E 2n = Ri + ... + R 2n = 2nE 1 + (z x + ... + z 2n )e 
2n - 1 1 2 , 



~ 2n{E l + 



(65) 



This result casts into the structure we guessed for the 
energy change of 2n pairs as induced by the Pauli exclu- 
sion principle between an even number N = 2n of similar 
pairs. 

(ii) If the number of pairs is odd, N = 2n + 1, we do 
the same but one of the Zi is now equal to zero at the 
lowest order in ^7. Eq. (61) is then replaced by 



~ ha-i 



1 



1 



1 



ai — a 2 a\ — 03 ai — a n 

1 ... + _L_ + J- (66) 



01 ai + a n 



a\ + a 2 2a\ 



So that the sum of a„'s just reads as Eq. (62) with n/2 
replaced by 3n/2. Consequently, we now find, instead of 
Eq. (64) 



1, I 



Y,^- IMn - l)/2 + 3n/2] 7/2/A 2 = n(n + ~hf 2 
i=i 1 



which ultimately leads to the expected result 
E 2n +i — Ri + ... + R 2n +i 



(67) 



2n I 2 

~(2n+l) [E 1 + — -| 



(68) 



F. Validity of the result 

Let us now come back to the validity of the procedure 
leading to Eqs. (65) and (68). 

(i) We first wish to note that these equations allow 
us to cast the energy of N Cooper pairs into the same 
general form given in Eq. (1), for odd and even TV, the 
interaction energy E int in Eq. (1) being such that 



IE. 



int 



111 

Poll 



Po 



-(l + 2t/ c ) 



1 

— + e v 
Po 



(69) 



1/po is the energy increase when adding one electron of 
each spin at the frozen Fermi sea level when the density 
of states is constant and equal to po- The number of 
pair states within the potential extension ft above | JFb) 
is Nq = poSl. The one-Cooper pair binding energy is e c 
while ty is the elementary contribution to the Cooper 
pair binding energy from each up and down spin empty 
state feeling the attractive potential. The result we find 
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leads us to write the energy of N electron pairs added to 
|F„) as 



E N ~ NE 1 + N ^ n l) E lnt 



N 



2 Uf 



N — 1 
2po 



1 



N- 



Nq 



(70) 



The N(N — l)/2 prefactor in front of -Ej„t clearly in- 
dicates that this part of the energy comes from interac- 
tion between two Cooper pairs. Due to the very sim- 
plified form of the reduced BCS potential V, interaction 
between Cooper pairs can only come from the Pauli ex- 
clusion principle between electrons. However, since the 
Pauli exclusion principle is TV-body by essence, higher 
order terms in N(N — l)(N — 2)... should appear in the 
iV-pair energy. A possible reason for their absence is the 
fact that we perform a 7 expansion to analytically solve 
Richardson's equations but we only keep the lowest or- 
der term. This leads us to think that the result given in 
Eq. (70) is valid provided that corrections induced by in- 
teractions are small compared to the main term, namely 
Ei = 2ef„ — e c . Since there are two physically different 
"main" terms, let us consider them separately. 

(i) As already highly remarkable, the first term of this 
iV-pair energy, 2 [Ne Fa + N(N - l)/2p ], is the exact ki- 
netic energy of iV free electrons with up and down spins 
when added at the Fermi level Sf 01 whatever N is - un- 
der the assumption of a constant density of states po- 
Indeed, electrons are fermions; the Pauli exclusion prin- 
ciple forces their energies to increase as Sf 01 £f q + 1/po; 

Sf + (N — l)/po; so that the total energy of N up 
and down spin electrons is equal to 



N-l 

^(normal) ~ X 

- 1 2^ 

p=0 



£F Q 



P_ 

Pa 



Ne Fn + 



N(N-l) 



2po 
(71) 

which exactly is the first term of Eq. (70). This shows 
that the procedure we have used to get the ./V-pair energy, 
with a validity a priori restricted to the dilute limit on 
the one-Cooper pair scale, still produces the exact free 
part of this energy, whatever N. Since this free part 
is nothing but the "normal" part of the superconductor 
energy, i.e., the one which does not cancel for V — 0, we 
can rewrite Eq. (70) as 



E N - E { ™ rmal) ~ —N I 1 — 



N-l 

Nq 



Eq.(72), a correction small compared to the main term 
imposes N <C Nq, i.e., N small compared to the total 
number of empty states from which the one-Cooper pair 
bound state is constructed. We note that this small N 
scale is not at all related to the one-Cooper pair scale 
which led to N < N c . 

(hi) From a mathematical point of view, the validity 
of Eq. (70) must be traced back to the validity of the 
7 expansion we have performed in Eq.(27) and that we 
have used, within one or two terms only, this 7 expansion 
being a priori valid for Zi small compared to 1. Due to 
Eqs. (63, 64), the sum of the z% real parts as well as 
the sum of the squares of the imaginary parts, increase 
as N 2 -f. Since 7 scales as 1/N C , as seen from Eq. (34), 
the fact that we use an expansion which is valid provided 
that all the |zj|'s are small, leads us to a priori restrict 
the validity of the obtained result to N small on the N c 
scale, or a power law of this scale, condition which is far 
more restrictive than TV <C Nq . Consequently, we do not 
a priori expect the result we have obtained in Eq. (70) 
to be valid in the dense BCS regime in which the number 
of added pairs is half the total number of pair states Nq 
in the potential layer £1. 

(iv) Nevertheless, the fact that the expansion we have 
performed, gives the correct normal part for the energy of 
N pairs, even when N ^> N c , is a strong indication that 
our result is most probably valid on a wider range. An- 
other strong support of this idea is the fact that this re- 
sult, when extended to a complete filling of the potential 
layer, i.e., to TV = Nq, gives an average pair binding en- 
ergy which reduces to zero. This is fully reasonable since 
all the available states are then occupied, so that noth- 
ing can change when the attractive potential is turned 
on since, due to its peculiar form, no exchange between 
pairs can exist. Consequently, the energy difference be- 
tween correlated and normal states must reduce to zero 
when N = Nq in agreement with Eq. (72). 

All this leads us to guess that, in spite of the mathe- 
matical procedure we have used to get the expression of 
En given in Eq. (70), which, to be valid, a priori restricts 
the number of pairs to iV <C N c , the validity range of our 
result is most probably far larger that this quite restric- 
tive condition which corresponds to an extremely dilute 
limit, the scale being the one-Cooper pair extension. 

Let us now compare the ./V-pair energy we have ob- 
tained within this canonical ensemble approach, with the 
standard BCS result obtained in the grand canonical en- 
semble. 



~-2Ne v [N n -(N-l)\ (72) 

where ^j"°™ ai ) [ s th e energy of these N pairs in a normal 
state, i.e., without interaction, the N added electrons 
having transformed the "frozen" Fermi sea \Fq) into the 
"normal" Fermi sea \F). 

(ii) If we now turn to the correlated part of the ./V-pair 
energy, i.e., the second term of Eq.(70) or the RHS of 



IV. LINK WITH THE BCS DENSE LIMIT 

A. The standard BCS result 

(i) One important characteristic feature of the BCS 
theory in its traditional form is that there is only one 
energy scale associated to pairing: the energy difference 
A between the ground state and the first excited state 
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linked to "pair-breaking" processes. Nevertheless, the 
gap A is often used to reach some physical understand- 
ing for the ground state itself. The energy difference 
between the ground state and the normal state having 
the same number of noninteracting electrons is written 
in textbooks as 



U, 



(73) 



where A reads in terms of the potential extension and 
amplitude as A = 2uj c e~ 1 / p ° v = fte" 1 ^ ^. 

By comparing the gap with the single pair binding en- 
ergy given in Eq.(9), we note a factor of 2 difference in 
the prefactors, but also in the exponent, the latter being 
crucial in the weak-coupling limit because it makes the 
gap far larger than the single pair binding energy. This 
immediately leads us to reconsider the widely accepted 
connection between " Cooper pair energy" and A because 
this binding energy can hardly increase with pair num- 
ber up to the gap value due to Pauli blocking among the 
available states to form paired states. 

(ii) The usual way to understand the energy difference 
in Eq. (73) is to say that each pair which lies in a narrow 
shell A at the noninteracting electron Fermi level, has an 
energy of the order of A. The number of such pairs 



N A ~ Po A 



(74) 



is very small compared to the total number of pairs feel- 
ing the potential. This gives an energy difference be- 
tween superconducting and normal states of the order of 
A^aA = poA 2 . This understanding indeed reproduces 
the result of Eq. (73), but within a numerical 1/2 pref- 
actor only, which is not that easy to bring in. 

(iii) By noting that the energy difference between su- 
perconducting and normal states also reads 



U s - U n 



1 



:Po 



Po- 



(75) 



there is another way to understand this result. Indeed 
the BCS configuration corresponds to add poQ/2 — Nq/2 
pairs to the frozen Fermi sea \F ), in order to fill half of 
the potential layer (see Fig. lb). Consequently, in order 
to recover the standard BCS result given in Eq. (75), 
the average binding energy of these Nq/2 pairs must be 
taken equal to Qe~ 2 / p ° v . This quantity which is just 
half the binding energy of a single Cooper pair e c , ex- 
actly corresponds to the average pair binding energy we 
have found in Eq. (72), when N is equal to Nq/2. Since 
the normal state energy U n is nothing but jr;^ ormal ) a p_ 
pearing in Eq. (70), we see that the BCS result given in 
Eq. (75) is identical to the energy of N Cooper pairs we 
have obtained in Eq. (72). Although the validity of our 
result seemed at first to be mathematically reduced to 
the very dilute limit N <C N c , due to the expansion (31) 



we have used to get it and the fact that we have only 
kept the first order term of this expansion, it is quite re- 
markable to note that its extension to the dense regime 
gives the BCS result exactly, i.e., with all its numerical 
prefactors. 



B. Dense limit with a non-symmetrical potential 

Such an agreement is even more remarkable when we 
consider a potential configuration which is not symmet- 
rical with respect to the normal Fermi level, i.e., for an 
arbitrary filling of the potential layer, not nessesarily 1/2 
as in the usual BCS configuration. 

(i) Following Eq. (2.33a) in Tinkham's textbook^, the 
energy gap obeys the equation 

-e F+ e Fo 



PqV 



sinh 



-1 ( £ F + 



sinh" 



V A 

l I £ F — £F 



(76) 



In the symmetrical configuration, £p is equal to £f + 
f2/2, so that the two terms in the above equation are 
equal. The value of A in the weak coupling limit (poV <C 
1) is then readily obtained as 



A ~ ne 



-1/poV 



(77) 



In a nonsymmetrical configuration, Eq. (76) cannot be 
solved exactly. However, if the Fermi level Ef is located 
deeply inside the layer, in order to have both ef +fi — £f 
and Ef — £f much larger than A, the inverse hyperbolic 
sines can be replaced by logarithms. Under these as- 
sumptions, Eq. (76) can be rewritten as 



1 



PoV 



In 



(£ Fa +£1-£f) (£F - £F ) 



which then gives 
A ~ 2 [(e Fo + n 



A 2 



£f) (£f - £F 



lV2 -,-1/poV 



(78) 



(79) 



The usual BCS expression (77) for the gap A readily fol- 
lows from this general result when the potential is sym- 
metrical, i.e., for ef — £f + fi/2. 

(ii) If we now use the expression of the gap given in 
Eq. (79) to calculate the energy of N pairs, we find 



Pa 
2 



(ef — £f 



(e f - £f„) V {£F - £F ) +A 2 

2 



+ ( e Fo +£L-£f) - ( £F a 

v /(e i , +O-e F ) 2 + A 2 



fl-e F ) 



(80) 
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For Ep + f2 — Sf and ep — ep both much larger than A, 
the Taylor expansion of the RHS of Eq. (80) shows that 
this energy difference takes the same form as the one for 
a symmetric potential given in Eq. (73), namely 



U s - U n 



1 



ill 



with A now given by Eq. (79) instead of Eq. (77). Sub- 
stitution of Eq. (79) into Eq. (81) then gives 



1 



which also reads as 



U.-U n 

- [po (SF ~ £F )] 

sf — £f 



n 



u s -u n 

-2Ne v (N n - N) 



(82) 



(83) 



The first factor in the RHS of Eq. (82) is the number N 
of pairs in the potential layer. The second factor of Eq. 
(82) is the binding energy e c of one isolated Cooper pair 
while the last factor, which reduces to 1/2 for a symmet- 
rical potential, in agreement with Eq. (75), just corre- 
sponds to the " shrinkage" of the average binding energy 
resulting from the Pauli exclusion principle. Eq.(82) thus 
shows that, in the nonsymmetrical configuration too, it 
is quite natural to take as Cooper pair number, the num- 
ber of pairs in the potential layer, their average binding 
energy being linearly reduced, due to the Pauli exclusion 
principle. Eq. (83) then shows that this binding energy 
is proportional to the number of empty states (Nq — N) 
remaining in the potential layer. 

We also wish to note that this BCS-like grand canon- 
ical treatment can be extended to the dilute regime of 
pairs provided that their number still is macroscopically 
large, i.e., their density is finite. In this case, one has 
to supplement the gap equation by an equation for the 
chemical potential. The resulting expression of U s — U n 
again agrees with Eq. (82). Conceptually, this BCS anal- 
ysis is similar to the Eagles's models for the density- 
induced BEC-BCS crossover (see also Refj22r— ). Al- 
though the BEC-BCS crossover problem is currently con- 
sidered mostly within cold gase s 26 ' 27 , there are some in- 
dications that a density-induced crossover also has some 
relevance for correlated pairs in high-T c cuprate o 28 ' 29 . 



C. Cooper pairs as composite bosons 

From a general point of view, the energy of N inter- 
acting pairs can always be written as 



NE 1 + N{N n l) E mt (n) 



(84) 



where n = N/L D is the pair density, LP being the sample 
volume. The result of Eq. (70), obtained in the dilute 
limit, gives Ei nt {n — > 0) = Ei nt . Due to the reduced 
potential taken in the BCS theory, with a constant inter- 
action between electrons having opposite spin and oppo- 
site momenta, it is hard to physically expect higher order 
density contributions in the interaction term coming from 
the potential itself. A possible N dependence of Ei nt (n) 
is more likely to come from the Pauli exclusion principle 
which is JV-body by essence. Actually, the BCS result 
clearly indicates that this additionnal N dependence does 
not exist. It can be of interest to note that we have re- 
cently found a similar result for the ground state energy 
of N Frenkel excitons: terms proportional to third and 
higher powers in N are absent in the energy of N Frenkel 
excitons^. The Hamiltonian for Frenkel excitons^ has 
similarities with the BCS Hamiltonian because it only 
couples electrons to holes localized on the same atomic 
sites: in the same way, the interaction potential used in 
the BCS theory only couples up and down spin electrons 
having opposite momenta: in both cases, this is a one-to- 
one coupling between different fermions, so that Cooper 
pairs like Frenkel excitons are composite bosons with one 
degree of freedom only, in contrast to Wannier excitons 
which have two: the electron momentum and the hole 
momentum. 



As a last comment, we wish to mention that the inter- 
action term Ei nt has two contributions which have com- 
pletely different origins although both of them end with 
the same N(N — 1) dependence. These two different ori- 
gins nicely enlighten the fact that Cooper pairs are boson- 
like particles made of two fermions. Indeed, one part of 
Ei nt is associated to a kinetic energy change coming from 
electrons added to the frozen Fermi sea \Fq). The sum of 
their kinetic energy increases l/po+2/p +... + (N—l)/p , 
equal to N(N — l)/2po for each spin, just gives the nor- 
mal part of Ei nt . The N(N—1) prefactor of this free part 
comes from a sum of N different energies with a arith- 
metical progression. In contrast, the interaction energy 
also has a N(N — 1) contribution which is fully linked 
to the bosonic nature of the Cooper pairs, these pairs 
all having the same binding energy 2[Nq — (N — l)]ey 
coming from the number of empty states in the potential 
layer: the resulting energy change then is 2(N — l)ey = 
e c (N — 1)/Nq. When multiplied by N, this also leads 
to a N(N — 1) prefactor in the energy of N pairs, but 
this prefactor results from the sum of N identical en- 
ergy changes. This understanding nicely reveals the 
two aspects of Cooper pairs: The bosonic nature of the 
pairs shows up in the correlated part of Ei n t , while their 
fermionic nature appears in the normal part of the en- 
ergy. This makes Cooper pairs quite interesting com- 
posite bosons, definitely, worth to be tackled along the 
ideas we have recently developed in the composite boson 
many-body theory 2 ^, originally constructed for Wannier 
excitons. 
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V. CONCLUSIONS 

We have extended the well-known Cooper's model be- 
yond one pair and revealed the deep link which exists 
between this model and BCS superconductivity. While 
Cooper has considered one electron pair added to a frozen 
Fermi sea, in a region where a small attracting potential 
acts, we here consider a large number of pairs, this num- 
ber being increased one by one, starting from N = 1. Wc 
calculate the ground state energy of these N pairs result- 
ing from both, the standard BCS potential and the Pauli 
exclusion principle. This is done by solving Richardson's 
equations for N energy-like quantities, their sum giving 
the exact TV-pair energy. We determine these N quanti- 
ties analytically through an expansion in the dimension- 
less parameter associated to the inverse of sample volume 
which turns out to also be the inverse of the number of 
pairs from which single Cooper pairs would start to over- 
lap. By keeping the first order term in this dilute limit, 
we find that the energy of N pairs reads as the energy 
of N individual Cooper pairs, within a N(N — 1) correc- 
tion which makes the average pair energy varying linearly 
with pair number. The average pair binding energy is 
found to be proportional to the number of empty states 
feeling the potential. Although this result only uses the 
first order term of an expansion which is mathematically 
valid in the dilute limit on the one-Cooper pair scale, 
quite remarkably, this linear behavior stays valid up to 
the dense BCS regime, even when the attractive potential 
between electrons is not symmetric with respect to the 
normal Fermi level, or when the added pairs completely 
fill the layer where the potential acts. The physics of 
this interaction term can be traced back to the Pauli ex- 
clusion principle, the number of available empty states 
within the potential extension decreasing linearly when 
the pair number increases. Since this argument is not af- 
fected when going from the dilute limit to the dense limit 
on the one-Cooper pair scale, it may appear as physically 
reasonable to find that the energy change obtained for N 
small on this scale, extends to the dense BCS limit in 
which a large fraction of the states feeling the potential 
is occupied. 

The nontrivial result we here obtain, brings some new 
light on the concept of " Cooper pair number" in the stan- 



dard BCS result for condensation energy. The usual un- 
derstanding is based on the idea that the Cooper pair 
binding energy is of the order of the excitation gap A; 
these pairs are concentrated in a very thin shell around 
the normal Fermi level, its width of the order of A be- 
ing much smaller than the width of the layer where the 
attractive potential between up and down spin electrons 
acts. However, in this understanding, the "number of 
Cooper pairs" is more a qualitative concept than a rigor- 
ously defined quantity. The understanding we here pro- 
pose uses the very basic fact that correlations exist be- 
tween all the electrons with opposite spins and opposite 
momenta lying in the potential layer. The binding energy 
we find for such pairs in the TV-pair configuration is much 
smaller than A, being half the energy of one isolated pair 
in the case of a BCS potential extending symmetrically 
on both sides of the normal Fermi level. This binding 
energy reduction entirely comes from the Pauli exclusion 
principle which reduces the number of states available to 
form bound pairs. This makes the average pair binding 
energy ultimately proportional to the number of empty 
states. Such a link between an isolated Cooper pair and 
the iV-pair configuration is going to be of importance for 
considering superconductivity in the framework of the 
many-body theory for composite bosons we have recently 
developed, Cooper pairs definitely being quite interest- 
ing composite bosons with characteristics rather differ- 
ent from the ones of Wannicr excitons which have up to 
now been the main application of this new many-body 
theory^. 
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